Boosted three-dimensional black-hole evolutions with singularity excision 
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Binary black hole interactions provide potentially the strongest source of gravitational radiation for 
detectors currently under development. We present some results from the Binary Black Hole Grand 
Challenge Alliance three-dimensional Cauchy evolution module. These constitute essential steps 
towards modeling such interactions and predicting gravitational radiation waveforms. We report 
on single black hole evolutions and the first successful demonstration of a black hole moving freely 
through a three-dimensional computational grid via a Cauchy evolution: a hole moving ~ 6M at 
0.1c during a total evolution of duration ~ QQM. 
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The accurate computational modeling of black-hole in- 
teractions is essential to the confident detection of as- 
trophysical gravitational radiation by future space-based 
detectors such as LISA and by the LIGO/VIRGO/GEO 
complex of ground-based detectors currently under con- 
struction. The sensitivity of these detectors will be sig- 
nificantly enhanced if accurate computer simulations of 
black-hole mergers can produce predictions of radiation 
waveforms 0. The Binary Black Hole Grand Challenge 
Alliance || was funded in September 1993 to develop 
the computational infrastructure necessary accurately to 
simulate the coalescence of black-hole binaries. The pri- 
mary objective of the resulting code will be the produc- 
tion of waveforms from binary black hole mergers. In this 
Letter we report on an important step towards achieving 
such simulations. 

A key difficulty in evolving black-hole spacetimes is 
handling the curvature singularity contained within each 
hole. The only viable means of accomplishing this over 
time scales required for binary coalescence appears to be 
black-hole excision: exclude all or part of the black-hole 
interior (and the singularity) from the computational do- 
main and evolve only the exterior region P,H. This 



is possible because, by definition, the region inside the 
hole cannot causally affect the exterior evolution. Black- 
hole excision has been implemented successfully in spher- 
ical symmetry by Seidel and Suen §], Scheel et al. §|), 
Marsa and Choptuik S, Anninos et al. |1, and Gomez 
et al. ju^; and in three dimensions (3D) by Daues Jll]] 
and by Gomez et al. [fl2|| . 

We are developing a general algorithm to perform a 
Cauchy evolution of Einstein's equations on a domain 
containing black holes with excised interiors. Prior sim- 
ulations of black-hole binaries have used coordinates in 
which the holes remain at fixed coordinate locations; for 
complicated motions, this will lead to undesirably con- 
torted coordinates. In contrast, our more general ap- 
proach allows black holes with excised interiors to move 
freely through the computational grid. Achieving this 
goal requires the synthesis of two key ingredients: a 
numerical scheme capable of stably evolving Einstein's 
equations on a domain containing moving boundaries 
(excised regions), and a set of gauge conditions that en- 
sure that coordinates remain well-behaved as black holes 
move through the grid. In this Letter we present a cru- 
cial step towards achieving a general black-hole evolu- 
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tion scheme: the first successful Cauchy computation of 
a single black hole freely moving through a 3D numerical 
grid. A different approach within the Alliance that uses 
a characteristic formulation |ljj has recently shown 3D 
black-hole evolutions to 1400M and success in moving a 
black hole less than one diameter. 

Black-hole excision is based on the fundamental idea 
that the event horizon of a black hole is a natural causal 
boundary. Unfortunately, an event horizon cannot be lo- 
cated without knowing the full future evolution of the 
spacetime. However, a related structure known as an 
apparent horizon can usually be located on a spacelike 
time slice using only the information on that hypersur- 
face. An apparent horizon is a topologically spherical 
spacelike two-surface on which the expansion of the con- 
gruence of outgoing null rays orthogonal to the surface 
vanishes. If the characteristic curves of all fields being 
evolved lie on or within the light cone, then the apparent 
horizon can be used as the inner boundary for a Cauchy 
evolution. The causal nature of this boundary implies 
that no explicit boundary condition need be imposed on 
the evolved quantities. 

The gauge freedom (coordinate freedom) of general rel- 
ativity allows considerable latitude in choosing how a 
computational solution evolves in time. Although gauge 
considerations cannot influence physics, they do deter- 
mine how the coordinates and computational grid points 
used to describe the solution relate to physical locations 
in spacetime. A poor choice of spatial or temporal gauge 
can lead to coordinate pathologies that ruin a numerical 
simulation. For example, the proper distance between 
two adjacent computational grid points might approach 
zero or grow without bound. It is not fully understood 
what constitutes a good gauge choice. 

An attractive choice [[l3|]§ll for describing a single black 
hole is the ingoing Kerr-Schild form [TJ|] of the Kerr met- 
ric: 

ds 2 = g lxu dx' J -dx v = (n^ + 2Hl /Jt l u ) dx^dx", (1) 

where fi,v run from to 3, — diag(— 1, 1, 1, 1), 
is a null four-vector, and H{x a ) is a scalar function. In 
this gauge, the coordinates are closely related to the null 
structure of the spacetime. Furthermore, the solution 
is time independent (or has a trivial time-dependence 
for a moving black hole) , the spacelike hypersurfaces ex- 
tend smoothly through the horizon, and gradients near 
the horizon are smaller than in several other coordinate 
choices. 

Equation (|lj) is form-invariant under Lorentz transfor- 
mations, so it can be used to represent either a non- 
translating or a boosted Kerr black hole. The 3+1 de- 
composition of the spacetime metric leads to a metric §ij 
and extrinsic curvature Ky that we use as initial data 
and for comparisons at later times ]15| |. The lapse func- 
tion, 



a = l/y/l + 2Hll (2) 

and shift vector, 

ft - IHkk, (3) 

are analytic functions of space and time that we im- 
pose as gauge conditions. For the time- independent 
Schwarzschild spacetime, H = M/r (M is the mass), 
Ifi = (l,Xi/r), and the apparent horizon coincides with 
the event horizon. In fact, the apparent and event hori- 
zons coincide in the boosted case as well. 

Schemes for excising the interior of a black hole from 
the computational grid typically require a superluminal 
shift vector in some region of the computational domain, 
and must cope with the lack of an explicit boundary con- 
dition on the excision boundary. The Alliance Cauchy 
evolution module implements an evolution scheme that 
is designed to provide a stable evolution for any choice of 
shift vector. A typical evolution equation has the form 

(d t -£p)T = ..., (4) 

where £p is the Lie derivative along the shift vector and 
T is an arbitrary spatial tensor. We can rewrite this as 

d Q T - (£f, - /m) T = ... , (5) 

where da = dt — /3 l di is a time derivative in the direc- 
tion normal to the spatial hypersurface and which, by 
definition, is centered in the light cone. This scheme re- 
moves from £p the advective term that potentially leads 
to evolution along a non-timelike direction. This guar- 
antees that a numerical evolution will be stable against 
instabilities produced by superluminal shifts. Note that 
the tensor components being evolved remain in the co- 
ordinate basis djdx % that is Lie dragged along the f 
direction, but the computational grid points at which 
these components are defined do not remain at constant 
values of x % . Instead, the grid points remain at constant 
values of the spatial coordinates x l that are Lie dragged 
along the unit normal to the spatial hypersurface. The 
evolved quantities are determined at the desired spatial 
coordinate locations x l by interpolation. This interpo- 
lation requires that the x % coordinate values be evolved 
along the t* 1 direction for each point where the evolved 
quantities are to be evaluated on a given time slice. It 
also obviates the need for a boundary condition at the 
excision surface. The algorithm can easily accommodate 
black holes that move through the grid. In this case, un- 
used grid points that had been in the interior of a hole 
can move into the exterior. These new points will be 
filled by the interpolation phase and will always be filled 
from data that is in the future domain of dependence 
of the previous time slice. Similar schemes were devel- 
oped by Seidel and Suen || , Alcubierre and Schutz ]lj| , 
and Daues (Till • Additional details on how this evolution 
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scheme is implemented can be found in Scheel et al. |7]] 
and in a forthcoming paper [ p~7[ that will describe the 
details of the Alliance Cauchy evolution module. The 
key point is that by splitting the evolution along into 
an evolution along the unit normal followed by an in- 
terpolation, the algorithm guarantees that the numerical 
evolution step is always taken in the center of the physical 
light cone. 

The Cauchy evolution module developed by the Al- 
liance is based on the standard (3+1), or ADM, decom- 
position of Einstein's field equations p8| , |l9| |. The module 
uses Cartesian coordinates and evolves the three-metric 
gij, and extrinsic curvature 2£y, on a 3D rectangular 
grid. The evolution equations are solved using an it- 
erative Crank-Nicholson differencing scheme. The con- 
straint equations are not imposed on the solution but 
are utilized as a diagnostic. 

We present results from two sets of black-hole ex- 
periments: evolutions of non-translating and of boosted 
Schwarzschild black holes. In both sets of experiments 
we work with a 3D nxnxm computational grid, and we 
impose Dirichlet boundary conditions on the grid faces 
throughout the evolutions. For the non-translating cases, 
this means freezing the outer boundaries to initial data. 
For the boosted case, this means resetting the exact solu- 
tion at the faces for each time step. At each time during 
the course of our evolution we impose the exact gauge 
conditions (^) and (||) as a function of four spacetime co- 
ordinates. We allow for a "buffer region" of p grid zones 
between the apparent horizon (where the inner boundary 
might have been placed) and the chosen inner boundary 
(excised region). 

For the case of a boosted black hole we discuss two 
runs; one with p — 3 buffer zones and the other with 
p = 5 buffer zones between the horizon and the inner 
boundary. Both runs use n = 33 and m = 65 with a 
domain of — (8/3)M to 8M in the z-direction and a do- 
main of — (8/3)M to (8/3)M in the x and y directions. 
The black hole is initially located at the origin and has 
a velocity v = 0.1c in the z direction. The evolution was 
carried out with a Courant factor of 1/4. Figure [j] shows 
a spacetime picture of g zz along the z-axis vs time for 
the p = 5 run. In this case we evolved the black hole 
for 61M, moving three hole radii. The evolution termi- 
nated when the leading edge of the black hole was within 
five gridpoints of the z = 8M face. At this point the 
differencing algorithms failed because of an insufficient 
number of points between the inner and outer bound- 
aries. In the p — 3 run the black hole evolved for 54M 
with coordinate stretching occurring at the trailing edge. 
We find that the addition of buffer zones enhanced the 
run and reduced coordinate stretching for the duration 
of the p — 5 run. This is consistent with the behavior 
of the non-translating runs described below. As the hole 
moved in each of the evolutions, many coordinate points 
emerged stably and smoothly from the excised region into 
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FIG. 1. Metric component g zz along the z-axis as a func- 
tion of time. The flat region that moves diagonally to the 
right represents the excised region (inside the black hole). 
Note that points at the trailing edge (left side) are smoothly 
updated as the hole moves towards positive z. Coordinate 
effects are seen to appear near the inner boundary. 

the computational domain. Figure ^ shows a spacetime 
picture of the normalized Hamiltonian constraint diag- 
nostic along the z-axis for the p = 5 run. Note that 
the Hamiltonian constraint is well behaved in the regions 
where grid points have emerged from the black hole. 

For the case of a non-translating black hole, we chose 
n = m = 49. The outer boundaries are placed at —AM 
and 4M in each direction. With a Courant factor of 1/4 
we can evolve a single black hole to 95M. This is an 
encouraging achievement comparable to the longest 3D 
single black hole Cauchy simulation (> lOOAf ; Daues 
0)- 

As in the moving hole case, the length of the evolu- 
tions is dependent on the placement of the inner bound- 
ary. We found that by setting the number of buffer zones 
to p = 0, 1, 5, and 9, we can run a non-moving hole to 
a maximum time of t max = 16M, 20M, 95M, and 82M. 
In all cases, the evolution terminates because the itera- 
tive Crank-Nicholson evolution scheme fails to converge 
to a solution. While the evolved solutions all deviated 
from the analytic solution, the behavior in each case was 
somewhat different and a definitive explanation of why 
each evolution fails requires further testing. The most 
likely causes of the late time problems are either coordi- 
nate effects and/or numerical instabilities due to choices 
of finite difference operators. One such coordinate effect 
can come from fixing the gauge via the analytic func- 
tions @ and (§). As numerical evolutions progress, the 
evolved data drift from the exact solution. When this 
happens, the choice of lapse and shift being used will no 
longer enforce the desired, underlying Kerr-Schild gauge 
condition. In this case, effects such as grid stretching can 
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straint, N = (R + K 2 - JT^yVQi?] + |A' 2 | + \K* 3 Ka\ r . 
along the z-axis as a function of time. The flat region shows 
the location of the excised part of the domain within the hole. 



occur and the coordinate location of the apparent hori- 
zon will no longer coincide with the analytic solution. For 
p = or 1, the evolution progresses rather smoothly and 
the computed location of the apparent horizon gradually 
moves inward. Eventually the apparent horizon passes 
within the computational inner boundary, so that the in- 
ner boundary becomes timclikc and therefore unsuitable 
as a boundary for black-hole excision. The evolutions 
terminate soon after this occurs. For p — 5 or 9, the 
evolution is more complicated. The relatively smooth 
growth in error with the apparent horizon moving in- 
ward no longer occurs (or occurs on a much longer time 
scale). In these cases, the dominant errors appear to be 
noise introduced at the Dirichlet outer boundaries. These 
errors propagate across the computational grid and ap- 
pear to be amplified in some way. It is possible that the 
analytic gauge conditions are causing this amplification; 
some unusual coupling with the inner boundary may also 
be responsible. Eventually, the geometry near the inner 
boundary becomes quite non-spherical. While we are no 
longer certain of the location of the apparent horizon at 
this point, we believe that the inner boundary again be- 
comes timelike prior to the evolution ending. 

The fact that the evolutions depend strongly on the 
placement of the inner boundary agrees with prior 3D 
work by Daues p"l[ ] but requires further study. Using a 
larger number of buffer zones forces the inner boundary 
deeper into the black-hole interior where an increasingly 
larger propagation speed is required for information to 
escape from the black hole. Thus the strong dependence 
of the evolutions on p could indicate some non-physical 
gauge-dependent quantities or numerical errors propa- 
gating faster than light. Resolution of this question must 
await at least a consistent set of convergence tests to un- 
derstand the effect of p as the discretization is refined. 
In spite of these concerns, the evolution to 95M is an 



encouraging result and we believe that improved gauge 
choices that utilize information from the evolution (e.g. 
the computed apparent horizon location, cf. Refs. f^,[il"||) 
will allow for much longer evolutions. 

In the evolutions of non-translating holes, each case 
above results in an unstable evolution. We believe that 
this is most likely due to coordinate effects. This be- 
lief is supported by evolving a region adjacent to a non- 
translating black hole. We consider a domain extending 
from —(3/2) M to (3/2)M in all directions and displace 
the center of the black hole along one axis by 2M. By 
this design the horizon passes through the center of the 
computational cube. Using resolution n = m = 49 and 
Dirichlet boundary conditions on all boundaries, we find 
that the solution settles down to a steady state because 
of numerical dissipation and perhaps because wave dis- 
turbances fall within the horizon. At low enough reso- 
lution, the errors in the analytic gauge conditions and 
Dirichlet boundaries were too large to allow a steady- 
state solution. As a rule, higher resolution cannot force 
a fundamentally unstable numerical scheme to be "more 
stable". We thus find this test to be highly suggestive 
that our interior evolution scheme is fundamentally sta- 
ble and our principal problems are with coordinate and 
boundary effects. 

All of the above results are very encouraging for the 
ultimate goal of evolving multiple black hole spacetimes. 
In both the boosted and non-translating evolutions, im- 
provements to the coordinate conditions are clearly nec- 
essary. Non-analytical coordinate conditions need to be 
developed for the boosted case and we are hopeful that 
it will be possible to formulate a general Kerr-Schild- 
like gauge condition that does not require the solution 
of elliptic equations. Other combinations of spatial and 
temporal gauge conditions may also work or may even 
be necessary. There is still considerable work to be done 
on the computational infrastructure required to support 
general binary black hole simulations, but it is clear that 
one of the keys to successful evolutions is understand- 
ing which gauge conditions are appropriate for multiple 
black holes. 
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